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Abstract. In a cellular mobile radio system an SDMA 1 component can be implemented for the reuse of radio 
channels physically incorporated by time, frequency or code slots. Since SDMA is based on the spatial separation 
of different users operating in the same channel, a DOA 2 sensitive channel allocation scheme is essential for max- 
imizing system capacity. In this paper we present the Eigenvector method, a computationally efficient algorithm to 
do this job. We also present simulation results on the Eigenvector method operating in a typical urban mobile radio 
cell. The simulated blocking probabilities are then used to predict the capacity increase which can be expected 
after adding an SDMA component to a conventional mobile radio system. 
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1. Introduction 

The basic idea of SDMA is the reuse of a radio channel - incorporated by an FDMA, TDMA 
or CDMA 3 slot 0- by K > 1 different users in the same mobile radio cell. 

On the uplink, the spatial separation of K signals can be done by exploiting the information 
supplied through an M > K element antenna array at the base. On the premise of roughly 
synchronized user bursts, the data sampled at the array can be used to jointly estimate the 
fast fading channel impulse responses relevant for all K user-specific uplink channels. These 
channel estimates can then be fed to a linear [1] [2] or a non-linear [3] [4] data detector 
yielding estimates for the symbols transmitted by each user. 

Alternatively, DOA based spatial filtering techniques can be employed to separate the K 
signals transmitted on the uplink. K single user data detectors operating on the separated 
signals can then produce estimates of the data symbols transmitted by each user [5]. 

On the downlink, the spatial separation of the users is more difficult, since the mobiles 
cannot necessarily be expected to do joint detection and adaptive interference cancellation: 

• The mobile terminals must be cheap, small and light. Preferably, the mobile station will 
not be equipped with an antenna array and sophisticated joint detection signal processing 
hardware. 

• In general, the propagation delay of the mobile radio channel varies from user to user. On 
the uplink, user-specific transmission timing can easily compensate for these delays. On 
the downlink, the K signals transmitted by the base station cannot be timed in a way that 
each user receives his signal synchronized with the interference from the signals bound 
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for the K — 1 co-users receiving in the same channel. Therefore, the mobile has little 
chance to effectively execute interference cancellation. 

For these reasons, the appropriate way to do spatial separation on the downlink is to prevent 
interference for all K users by beamforming at the base station antenna array. The beam- 
forming weights are supposed to be controlled by the spatial channel parameters estimated by 
means of the K user signals received on the uplink. 

In many mobile radio systems the channel impulse responses estimated on the uplink can- 
not be directly reused as beamformer inputs due to the frequency and/or time gap between the 
uplink channel and the downlink channel [6] [7] [8]. In fact, channel reciprocity is only given 
for the DOAs and the corresponding medium term averaged path attenuations. Therefore, the 
adaptive control of the beamforming weights has to be based on the uplink channel estimates 
averaged over the fast fading. 

The medium term downlink channel of each user k = 1 • • K can be efficiently described 
by means of the M x M spatial covariance matrix 

C k = T,%Al g a kq a» g . (1) 

The number of propagation paths between the location of the base station and the location 
of the user k is denoted by Q k . Each path q = 1 • • Q k is described by its average squared 
transmission factor A kq and the array steering vector a kq incorporating its DOA, represented 
by an azimuth-elevation pair (ir kq , 0 kq ). Efficient algorithms named Unitary ESPRIT and 2D 
Unitary ESPRIT to estimate these parameters in real time were presented in [9] and [10]. 

The spatial covariance matrices C\ • • Ck do not contain any information about the fast 
fading relevant for downlink transmission. Therefore, the average SNIR 4 necessary for each 
downlink receiver will be higher than the average SNIR the base station antenna array can 
cope with during uplink reception. 

This implies that in contrast to most conventional mobile radio systems, the SDMA down- 
link, not the SDMA uplink, will be the critical link. Hence, the SDMA capacity increase 
gained by reusing the resources within a cell has to be evaluated according to the average 
number of users which can spatially separated by beamforming on the SDMA downlink. 



2. Spatial Separability 

On the downlink, the spatial separation of K different users can only be managed in a robust 
way, if the DOAs of all users in one channel are not too close to each other. Otherwise, 
downlink transmission will face severe problems, which can be illustrated by the following 
example: K = 2 users, each characterized by a single propagation path with identical atten- 
uations, have to be separated by means of an M —4 element uniform linear antenna array 
(ULA). In case (a) the azimuths , \j/ 2 of the two corresponding DOAs are quite different 
from each other (V^i = 4-30° , ty 2 = -30°), whereas in case (b) they are almost identical 
=+3°, V 2 = -3°). 

The beampatterns created by two optimized weight vectors w\ and w 2 yielding an average 
SNIR of 10 dB for both users are depicted in Figure 1. In case (b) the steep slopes of the 
patterns at the DOAs of the users (±3°) make two problems obvious: 
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1. The performance of the downlink beamformer is extremely sensitive to DOA estimation 
errors. 

2. The DOAs of the users are far away from the maxima of the corresponding beampatterns. 
This waste of electromagnetic power will result in unnecessary CCI in neighboring cells. 

In this context, it makes sense to call the case (a) a "spatially well separable" scenario and the 
case (b) a "spatially badly separable" scenario. 

The problems occuring in scenarios like (b) can be avoided by combining SDMA with 
at least one different multiple access scheme like FDMA, TDMA or CDMA which can sup- 
ply the system with a number L of separate channels. A DOA sensitive channel allocation 
algorithm can then assign spatially badly separable users to different channels. 

3. DOA Sensitive Downlink Channel Allocation 

Let us assume that prior to a new user's channel request there are K (l) users operating in each 
channel / = 1 • • • L. The numbers K (l) are not necessarily equal. In order to avoid trivial 
solutions there are no vacant channels assumed (K (l) =/ 0). Altogether J — 1 = A" (,) + ■ - * 4- 
K {L) users are on air. 

After adding a new user to the system, L J new user-channel combinations will be possible, 
provided there are no restrictions concerning reallocations of the J — 1 active users. In this 
case, the number of combinations to be checked will be prohibitively high even in small 
systems (e.g. 7 • 10 !1 for L = 7 and K — 14). Another reason, why we will not consider any 
user reallocations during the allocation procedure, is the additional signalling traffic caused 
by intracell handovers. 

Therefore, we will solve the SDMA channel allocation problem by a two-step procedure: 
First, finding the channel / opt e { 1 • • • L } with best spatial separability for the new user; 
second, evaluating the spatial separability in that channel and then decide whether the new 
user will be allocated to the chosen channel or whether his channel request will be rejected. 

The quality of any SDMA channel allocation scheme has to be evaluated according to how 
far the following goals can be achieved: 

1 . Maximizing system capacity by maximizing the average number of users which can be 
accommodated in L channels. 

2. Ensure high probability for robust downlink communication by maximizing robustness of 
the beamformer against parameter estimation errors. 

3. Minimizing CCI in neighbouring cells by minimizing the average RF power emitted by 
the base. . 

A computationally efficient algorithm doing this job is given by the Eigenvector Method 
(named "Quick SB algorithm" in [11]). It is based on the considerations presented in the 
following sections. 

3.1. The SDMA Downlink Beamforming Problem 

Let us assume we want to supply one user, characterized by the receiver noise N (composed 
of thermal receiver noise and CCI from neighboring cells) and the spatial covariance matrix 
C, with a given signal-to-noise-and-interference ratio SNIR by applying the complex weight 
vector w at the M element base station antenna array. The weight vector w will be chosen in a 
way that it minimizes the downlink transmit power P that has to be emitted at the base. 
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Figure 1. Beampatterns created by a 4-element ULA for a spatially well separable scenario (a) and a spatially 
badly separable scenario (b). 
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The downlink transmit power P is proportional to the squared length w H w of the weight 
vector applied at the array, whereas the RF receive power S at the mobile antenna is pro- 
portional to the term w H Cw. Therefore, the beamforming problem to calculate the optimum 
weight vector w can be mathematically put as the following constraint optimization problem: 

minimize r p = / i subject tQ s = w H Cw = N ■ SNIR. (2) 
w 

The solution n> of the above problem is proportional to the dominant eigenvector «(0 of the 
spatial covariance matrix C, i.e. the eigenvector corresponding to the largest eigenvalue X{C) 
of C. Hence, the minimum downlink transmit power is given by the product X~ X (C)-N - SNIR. 

Generalizing this result to the case of K users (indexed by (•), - • ( ) K ) being accommo- 
dated in K separate channels leads to the minimum RF transmit power 



* = 1 



K N k SNIR 



Let us now consider the SDMA case with K users operating in the same downlink channel. 
The beampatterns to separate the users k = 1 • • K from each other on the downlink will be 
produced by the weights w\- — wk. Generalizing (2), the corresponding SDMA beamforming 
problem turns out to be an optimization problem with K constraints which can be put as 
follows: 



minimize 



{p = f>f*J subjectto 5f^ = A fc + ^^C*Hr/ V *=!...*. (4) 
1 J £1 



3.2. Measuring the Spatial Separability 

If SNIR, Nu--,N K > 0 holds and all spatial covariance matrices C k are non-negative def- 
inite (following from the definition (1)), the downlink transmit power P solving the SDMA 
beamforming problem (4) will never be lower than the minimum power P min defined in (3). 

The case P = P mm is characterized by the dominant eigenvector of each user's spatial 
covariance matrix being orthogonal to all other users' spatial covariance matrices: 

C k S(Q) = 0 V k = 1 . . . K , b=/ k = 1 . . . K. (5) 

Plugging these conditions into (4) shows that in this case all weight vectors w k will be pro- 
portional to the dominant eigenvectors w(Q), like in the single user case (2). The downlink 
beamformer only produces useful signal power but no cross interference, since ytf C k w\ = 0 
holds for all h=/ I. Therefore, it makes sense to refer to the case incorporated in (5) as "perfect 
spatial separability". 

Simulations have shown that for "imperfect spatial separability" there is a connection 
between the angular distance of the users and the difference between the downlink transmit 
power P and the theoretical minimum P^ xn . As an example, let us assume K = 2 users have 
to be separated on the SDMA downlink by beamforming with an M — 8 element ULA. The 
channels of both users can be described by means of a single DOA each. The corresponding 
attenuations (p x = p 2 = 0 dB) and the elevation angles (0 X = 0 2 = 0°) are identical. Both 
DOAs are considered to be made up of very large numbers Q\ and Q 2 of propagation pathes 
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with Gaussian distributed azimuths. The mean azimuths and \fr 2 are different, whereas the 
corresponding standard deviations ("angular spreads" [12]) are identical (s\ = s 2 = s). 

Figure 2 shows the ratio P / P min , with the first user constantly looking at the array boresight 
(i/rj = 0°), and the azimuth j// 2 of the second user varying from —90° to -{-90°. Ignoring 
the sidelobes in Figure 2, the ratio P/P m i n tends to increase, if the DOAs of the two users 
move closer to each other. Finally, if the DOAs are just too close (or even identical), the 
problem (4) does not yield any solution at all, i.e. the users are spatially no longer separable 
by beamforming. 

Therefore, it makes sense to use the ratio r) = P/P mm , which from now on will be refered 
to as the "cross interference loss", as a measure for the spatial separability of a scenario. Doing 
this, rj = 0 dB means optimum spatial separability, whereas t] — ► oo refers to the case of no 
spatial separability at all. 

Note, that since spatial covariance matrices are the input parameters of the beamforming 
problem (4), the definition of the cross interference loss is not restricted to single DOA sce- 
narios or discrete DOA scenarios. The cross interference loss can rather be calculated for any 
narrowband scenario and any antenna array configuration. 

3.3. Calculating the Cross Interference Loss 

The cross interference loss r) depends on the downlink transmit power P and is, therefore, 
only known after the SDMA beamforming problem (4) has been solved. Unfortunately, this 
problem is numerically tricky, since both the target function and the constraints are non-linear 
functions of the variables w\ • • . 

Computationally inexpensive algorithms yielding approximations of the solution of the 
SDMA beamforming problem have been presented in [13] and [14]. The computational costs 
of these algorithms are basically caused by K generalized (Hermitian) M x M eigenvalue 
decompositions (GEVDs). If a DOA sensitive channel allocation scheme has to estimate the 

cross interference loss r\ in each channel / = 1 - - L, altogether L 4- K (l) H h K (L) GEVDs 

have to be executed for each user request. The resulting computational effort might be critical 
in real time, especially in large systems with a high number L of channels. 

Therefore, we will present a new method to quickly calculate an approximation of the 
solution of (4). We will refer to it as the "eigenvector approximation", since the basic idea 
of this algorithm is to restrict the weight of each user in a way that it is proportional to the 
dominant eigenvector of the corresponding spatial covariance matrix: 



This way the degrees of freedom are reduced to a great extent, since the complex- valued M x 1 
variables w\ • ■ w K will be replaced by the real-valued 1 x 1 variables Pi • • • Pk- 

Plugging the extra constraints (6) into (4) reveals that the originally non-linear constraint 
minimization problem melts down to an unconstraint linear problem: 



Solving (7) yields the variables P\ - Pk and the approximated downlink transmit power 
P = P, + - - • + P K . Like this, the computational costs for each user request are basically 
resulting from the need to compute L dominant eigenvectors (one in each channel / = 1 • • • L) 




(6) 




(7) 
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Figure 2. The cross interference loss r) versus the azimuth ^2 of the second user. 

for the matrices cf of the new user. Note that in mobile radio systems characterized by spatial 
covariance matrices that are almost identical for all channels (which is true for GSM and most 
of the other present systems), only one dominant eigenvector has to be calculated for each 
user request. 

Figure 3 shows the approximated cross interference loss fj = P/P m in for the scenario 
described in Section 3.2. Comparing Figures 3 and 2 exhibits that the eigenvector approxima- 
tion works particularly well for spatially well separable scenarios, characterized by a small 
cross interference loss r). This is not surprising, since in the case of perfect spatial separability 
(n — 0 dB) the eigenvector approximation yields the exact solution to the original SDMA 
beamforming problem (see Section 3.2. 

Evaluating Figure 3 for scenarios with high cross interference reveals that the eigenvector 
approximation produces more "pessimistic" results about the spatial separability of scenarios 
than an exact algorithm would do. But since in practice only spatially well separable sce- 
narios with small cross interference losses tj are relevant for downlink communication, the 
eigenvector approximation suggests to be an efficient tool to evaluate the spatial separability 
of a mobile radio scenario. 

3 A The Eigenvector Method 

Finally, the Eigenvector Method, a channel allocation scheme based on the considerations in 
the previous sections, is given by the following procedure: 
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Figure 3. The eigenvector approximation fj of the cross interference loss versus the azimuth \jf 2 of the second 
user. 



1 . Choose appropriate values for the desired signal-to-noise-and-interference ratio SNIR and 
for the maximum tolerable cross interference loss fj max . 

2. Estimate the spatial covariance matrix C of the new user requesting for a communication 
channel. 

3. Compute a dominant unit eigenvector m(C) of C. 

4. For all channels / = 1 - • - L do: 

• Estimate the receiver noise JV (/) at the location of the new user's mobile in the specific 
channel /. In general, this value N (l) predominantly results from CCI from neighboring 
cells and has to be measured by the mobile and communicated to the base station. 

• Assume the new user will be allocated to the channel /, so that the channel will have 
K — K il) + 1 users indexed by (-)i •••()*• With all spatial covariance matrices 
cfP • • • the corresponding dominant unit eigenvectors mCC^ 0 ) • • • "(C^) and the 
noise powers A^, (/) • • • known to the base station, solve the real-valued K x K system 




(8) 
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Figure 4. Traffic model for an SDMA system with L channels. 

• If the system does not have a (unique) solution or any entry is non-positive, the K 
users in the channel / will be considered spatially unseparable. Otherwise, calculate the 
eigenvector approximation rj il) of the cross interference loss: 

£ Pt l) 

rjM = -J=l . (9) 

SNIR Yi^'H^)- N ( k n 

5. Select the optimum channel / opt G { 1 • « • L } according to the smallest approximated cross 
interference loss. 

6. If ^ (/ °P l) is larger than the threshold ?7 max , then 

reject the new user (since the system is supposed to be fully loaded), 
otherwise 

allocate him to the channel / opt and calculate the corresponding 
beamforming weights: 

if = yf^F^). (10) 



The parameters SNIR and fj max must be both chosen according to the system parameters 
and the type of radio channel. We assume that for present mobile radio air interfaces and 
Rayleigh fading channels an SNIR between 10 and 20 dB and an fj max between 1 and 5 dB will 
guarantee robust downlink communication for every user in every channel. 

4. Capacity 

We will define the capacity of a mobile radio system in a traffic theory like manner: 

Capacity is the traffic A = XI an SDMA system can support without exceeeding a 
maximum blocking probability B during the channel allocation procedure. 
In this context the calling rate (in calls per second) is denoted by X, whereas T designates the 
average duration of a call (in seconds). 

The corresponding traffic model is depicted in Figure 4: Each state in the Markov chain is 
characterized by the number J of users in the system. The transition from a state J = j to 
a lower state J = j — 1 results from the termination of a call by one of the J users. Hence, 
the corresponding transition rate is given by the ratio j/ T of the number of active users to 
the average call duration. The transition from a state J — j to a higher state J = j + 1 is 
triggered by a user request resulting in a successful channel allocation. As long as there are 
still free channels available (j < L), the corresponding rate is identical to the calling rate 
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X. If all channels are occupied with at least one user (j > L), there is a chance the channel 
allocator rejects the new user, represented by the rejection rate r(j). 

For M > 1, there is no way of describing the rejection rates r(J) and tne resulting 
overall blocking probability B by means of analytic formulae. Therefore, we had to resort 
to simulations in realistic SDMA scenarios to estimate B as a function of the traffic A = XT, 

5. Simulation Results 

Our SDMA capacity predictions were based on a Monte Carlo simulation of a single mobile 
radio cell. The statistics of the parameters denning the radio channels between the users and 
the base station were chosen in compliance with the results of the 2D channel measurements 
carried out in the city of Munich in 1995 [15]. 

We are assuming a ring-shaped cell which has an SDMA base station in its center equipped 
with a uniform linear M element antenna array. The maximum distance from any user to 
the base station is given by the outer ring radius 5 km, the minimum distance by the inner 
ring radius 0.1 km. The user locations are uniformly distributed in the ring and independent 
from each other. Assuming a typical urban area, the average attenuation corresponding to the 
distance r k of the user k to the base station can be approximated by p k = 40 lg(r k / 5m) dB. 
The corresponding attenuations are log-normally distributed with the Suzuki parameter S = 6 
dB and the average p k . 

The propagation paths of each user k are assumed to consist of either one group (G k = 1, 
e.g. LOS 5 or one dominant reflector) or two groups (G k = 2, e.g. LOS plus one dominant 
reflector or two dominant reflectors). The probabilities for both cases G k = 1 and G k = 2 
are 0.5 each. Each group g = 1 • • • G k is assumed to consist of R kg e { 1 - • • 50 } propagation 
paths. 

The azimuths xjr kgr and elevations 0 kgr of all paths r = 1 ■ • ■ R kg are both Gauss-distributed . 
with the means ^r kg and 0 kg and the standard deviations = s$ — 5°. The mean azimuths 
lfr kg are uniformly distributed in the range [-180°; +180°], whereas the mean elevations are 
constant (0 kg = 0°). 

The noise powers N ( k n are assumed to be made up by CCI from log-normally shadowed 
base stations 1 5-25 km away from the users. This interference model corresponds to a hexag- 
onal non-SDMA cellular mobile radio network with a frequency reuse factor r — 4. The 
Eigenvector method operates with the signal-to-noise-and-interference ratio SNIR = 10 dB 
and the threshold ?7max = 

3dB. 

The dependence of the blocking probability B on the traffic A = XT is depicted in Figure 5, 
with a constant number L of channels (top plot) and a constant number M of antennas (bottom 
plot). 

Figure 6 shows the number L of separate channels an SDMA base station needs to support 
a given traffic A in the cell. The number of antennas was varied from M = 1 to M = 16. 
The number of calls simulated for each point in the plot is 10000. The tolerable blocking 
probability was 5 = 1%. 

Not surprisingly, Figure 6 can be interpreted in a way, that the higher the number M of 
antennas, the lower the number L of channels necessary to support a given traffic A. As an 
example, consider the traffic A = 30 erl: A conventional system (M = 1) like GSM needs 
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Figure 5. Blocking probability Z? versus traffic A. 
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Figure 6. Supportable traffic A versus the number L of channels with a maximum blocking probability B = 1%. 

L\ = 42 channels in order not to exceed the blocking probability B = 1%. An SDMA system 
with M = 8 antennas needs L 8 = 26 channels and one with M = 16 antennas needs L\ 6 = 14 
channels. The corresponding results for A = 60 erl are: L x — 73, L 8 = 32 and Lj 6 = 24. 

6. Conclusions 

In this paper, SDMA capacity is defined as the traffic A a system can support without exceed- 
ing a maximum blocking probability B during the channel allocation procedure. The capacity 
was estimated by simulating the performance of a specific SDMA channel allocation scheme, 
the Eigenvector Method^ in a realistic urban mobile radio cell. We assumed that channel 
allocation and, hence, system capacity is limited by the maximum number of users which 
can be accommodated on the downlink in a robust way. 

As a result, the simulations yielded the number L of channels (i.e. FDMA, TDMA or 
CDMA slots) a base station must provide to be capable of handling a given traffic A without 
exceeding the blocking probability B = 1% (see Figure 6). The capacity increase L X /L M over 
a conventional system with a single antenna is rather dependent on the number M of antennas 
than on the traffic A, as shown in Table 1 which has been extracted from the plot shown in 
Figure 6. 

The SDMA capacity increase L\/L M refers to the RWC gain from managing more than 
one user in each channel of a mobile radio cell. Note that employing SDMA features in all cells 
of a mobile radio network will result in an even higher increase in total system capacity. The 
reason is that downlink beamforming with antenna arrays in each cell produces less CCI in 
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Table 1. The SDMA capacity increase 
L\/Lm over a conventional system depend- 
ing on the traffic A and the number M of 
antennas. 



l x il m 


/4 = 20 


A =40 


A =60 


M = 1 


1.00 


1.00 


1.00 


M =4 


1.49 


1.66 


1.73 


M = 8 


2.12 


2.24 


2.38 


M = 12 


2.53 


2.63 


2.74 


A* = 16 


2.75 


2.95 


3.06 



neighboring cells than omnidirectional transmission with single antennas. Therefore, system 
capacity can be further increased by reducing the frequency reuse factor employed in the 
network. 
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